A New Class of Discretization Methods for the Solution of Linear Differential–algebraic Equations with Variable Coefficients

We discuss new discretization methods for linear differential–algebraic equations with variable coefficients. We introduce numerical methods to compute the local invariants of such differential–algebraic equations that were introduced by the authors in a previous paper. Using these quantitities we are able to determine numerically global invariances like the strangeness index, which generalizes the differentiation index for differential–algebraic equations that in particular include undetermined solution components. Based on these methods we then obtain regularization schemes, which allow us to employ general solution methods. The new methods are tested on a number of numerical examples.